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This article describes some recently obtained results on the low-energy properties of the "Kondo quartet" model 
of two spin- 1/2 impurities interacting with two channels (flavours) of conduction electrons. We shall particularly 
emphasize the connections between conformal field-theory methods and bosonisation approaches, which are first 
illustrated on the example of the single-impurity, two-channel Kondo problem. This article is dedicated to the 
memory of Claude Itzykson, and will appear in the Proceedings of the Conference "Advanced Quantum Field 
Theory", held in La Londe Les Maures in Sept, 1996 (Nucl. Phys. B, Proc. Supp., V.Rittenberg, J.Frolich and 
A.Schwimmer eds.). 



1. The two channel Kondo model 

Two methods leading to a low-energy solution 
of the single impurity, two channel Kondo model 
jl] have been recently put forward. One, due to 
Affleck and Ludwig 0,0] , uses the general frame- 
work of (boundary) conformal field theory (CFT). 
The other, due to Emery and Kivelson 0,^, relies 
on abelian bosonisation and establishes an exact 
mapping onto a resonant-level model which re- 
duces to free fermions for a special value of the 
Kondo coupling J^ . This is analogous to the 
Toulouse limit [g[ of the single-channel case. The 
aim of this section is to provide the reader with 
an introduction to the two channel Kondo physics 
by reviewing partly these methods and, more im- 
portantly, by establishing the precise connection 
between the two approaches. This will prove to 
be very useful in the study of the corresponding 
two impurity problem. 

1.1. Model and Symmetries 

We formulate the model in terms of left-moving 
chiral fermions ipia{x) on the full axis — oo < x < 
-|-oo, and of a spin-1/2 impurity spin S placed at 
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the origin, i = 1, 2 is a channel index, and a =|, | 
is a spin index. The hamiltonian is written as: 

H = ^VF E» /-r dx^lix)^^P,^ix) (1) 

+ JzS,JM + J{S^MO) + SyJyiO)) (2) 

In these formula, 

Jix) = ELiEa/3V'L(a^)^V'»/3(a;) denotes the 
total conduction electron spin current at position 
X. A Kondo coupling which is anisotropic in spin 
space has been considered, for reasons which will 
become clear below. The anisotropy J — Jz is 
actually irrelevant at the non-trivial fixed point. 
As demonstrated by Nozieres and Blandin m 
and widely discussed since, both the free fermion 
(J = 0) and the strong-coupling (J = +oo) fixed 
points are unstable under renormalization. The 
low-energy physics is controlled by an interme- 
diate coupling fixed point, which has non Fermi 
liquid properties Ximp ~ C/T ~ Inl/T. 

The interacting hamiltonian has, in the 
isotropic case Jz = J, a global symmetry 

SU{2)spin O SU{2)fiavour «> U{l)charge- This 

global symmetry group is smaller than the sym- 
metry of the non-interacting theory with Jz — 
J = 0, which has a t/(4) = 5(7(4) (g) C/(l) in- 
variance mixing all four species i, a. This invari- 
ance is partially broken by the Kondo interaction. 
In fact, as explained in detail in P|, the non- 



interacting theory has a full SO {8) global symme- 
try including transformations which mix flavour 
and charge, and the interacting theory itself has 
a larger global symmetry than the spin-flavour- 
charge one, namely an SU{2)spin ® SO{5)cf in- 
variance. We shall make use of this larger sym- 
metry later in this paper. 

The basic assumption of the CFT approach is 
that, at a fixed point, these global symmetries 
are promoted to local conformal symmetries. The 
symmetry algebra becomes the Kac-Moody alge- 
bra built on the global symmetry group of the 
interacting theory, namely: 



SU2{2)s(E)SU2{2)f(g)U{l) 



(3) 



Here, 51/^(2) stands for the level-fc SU{2) Kac- 
Moody algebra H. The Kac-Moody symmetry 
puts strong constraints on the theory. In partic- 
ular, it dictates the general form of the finite-size 
spectrum at any fixed point of the model as: 



LAE _ j(j + 1) 



J/(J/ + 1) 



nvp 






In this expression, AE = E — Egg are the ex- 
citation energies, vp is the Fermi velocity and 
L the radial length of the system. The quan- 
tum number j labels the representation of the 
5'C/2(2) j„ algebra and takes the possible values 
j = 0, 1/2, 1. Similarly, jf = 0, 1/2, 1 labels the 

representation of SU2{2) fi^^^^^ and Q is a (pos- 
itive or negative) integer corresponding to the 

charge of U{1)^. m,mf,mc = 0, 1, 2,3, ...'s are 
integers corresponding to towers of states built 
on these primary states. The degeneracies are 
also specified uniquely in terms of these quantum 
numbers. 

The general form ^ dictates the possible en- 
ergy levels, but does not specify which set of 
quantum numbers j, j/, Q actually appears in the 
spectrum at a given fixed point. The strategy 
followed by Affleck and Ludwig is first to identify 
the appropriate selection rules for the free fermion 
fixed point (J = 0). For this fixed point, all en- 
ergy levels can be built following the Pauli princi- 
ple. A major insight is then that the spectrum at 
the non-trivial infra-red stable fixed point can be 
obtained from a 'fusion principle'. Specifically, 



the spectrum is obtained by acting on the pri- 
mary operator in the spin sector associated with 
a given state of the free-fermion fixed point with 
the operator (j = 1/2) of the SU2{2) -^ algebra 
(leaving jf and Q unchanged). Each admissible 
values of j — jpp at the non-interacting fixed 
point then gives rise to at most two admissible 
values at the new fixed point, according to the 
fusion rule of CFT: 

j = 0^3' = 1/2 
j = l/2^/ = 0and/ = l 

J- = 1 ^ / - 1/2 (5) 

Physically, this fusion principle reflects the (par- 
tial) screening of the impurity spin by the conduc- 
tion electrons. The conformal towers appearing 
in the resulting spectrum are detailed in Ref. pH . 
This construction can be extended to find the op- 
erator content at the non-trivial fixed point (by a 
double fusion). 

1.2. Bosonisation 

The decomposition of a fermion ^ia{x) into 
spin-flavour-charge degrees of freedom can be 
made explicit in two possible ways. The first is to 
make use of 'non-abelian bosonisation' which has 
the advantage of respecting all global symmetries 
explicitly. This route is the natural one within 
the CFT approach. The fermion field is written 
as: 



In this expression, 17+, h^ 



(6) 

are Wess-Zumino- 
Witten fields belonging to the fundamental rep- 
resentation of the SU2{2)s and SU2{2)f spin 
and flavour algebra, respectively, and $c is a 
charge boson. Alternatively, one may take a more 
'pedestrian' (and often more explicit) route by us- 
ing abelian bosonisation, even though the global 
SU{2) invariances are no longer manifest. This is 
the starting point of the approach due to Emery 
and Kivelson Q , who introduce four boson fields 
such that: 

i,,^{x) = ^±=e-^*'"(-) (7) 



y/2TTao 



where: 



dx'Hiaix') - (l)iaix)} 



(8) 



with: 



[(l)iaix),Ilj,3{x')] = iSijda/sSix - x') 



(9) 



Our conventions for boson fields are such that 



an operator O 



o»9* 



has dimension A 



Oil 



Oil 



(t>f = o{(t>i^ +' 



Osf 



5l| 



»u 



'u 



'u 



'u 



02 T 



02 1 



02 1 



02 T 



»2iJ 



»2iJ 



»2iJ 



'2iJ 



(10) 

(11) 

(12) 
(13) 



g /2 {i.e its correlation function behaves as < 
O(0)O(x) >- 1/x^^ for |x| ^ +oo). 

Linear combinations are then formed, cor- 
responding to charge <j)c, spin cfig, flavor 0/ 
and spin-flavor (j)sf degrees of freedom (similarly 

1 

2 
1 

2 

In term of these fields, the components of the to- 
tal spin current read: 

•^^ ^ ;r^'^'''5'^s'^°**s/ , Jy = -^sin^scos^sf 

J - L§±^ (14) 

Only $s and <i>s/ enter these expressions. This is 
expected since these fields correspond to the sum 
and the difference of the bosons associated with 
spin degrees of freedom for each channel, namely 
$i-|- — $i| and $21 ~ <&2X- Furthermore, note that 
among the four independent real combinations of 
these fields, only three appear: sin^sf is not in- 
volved in the total spin current. That three real 
fields enter the current reflects the fact that the 
SU2{2) Kac-Moody algebra has central charge 
c = 3/2. These fields have dimension 1/2 and can 
alternatively be considered as three real fermionic 
fields (Majorana fermions) x" i o = a;, y, z: 

X^ = sin^s , X^ = cos<^s , X"" = cos<^sf (15) 
in terms of which the total spin current reads: 

J"" = ^eabcx'x' (16) 

This representation of the S'C/2(2) algebra in 
terms of three Majorana fermions is well known 
in CFT [0 , and will prove to be most useful below 
in the solution of the two impurity problem. 



The hamiltonian is easily written in terms of 
the boson fields. The free (kinetic energy) part 
reads: 



while the Kondo interaction takes the form: 



-.d-i-s 



^\S^cos<^s{0) + Sysin<i>siQ)]cos^sf{0) (19) 



1.3. Fixed-point boundary condition: the 
fusion rule hypothesis derived from 
the bosonisation approach 

In order to proceed with the analysis of Hint , 
Emery and Kivelson [Q observe that $s can be 
eliminated from the terms involving S^ and S'^, 
by performing a rotation of the fields along the 
S^ direction, i.e acting with the operator: 



Z.^ = exp(iS'^$s(0)) 

The transformed hamiltonian reads: 



(20) 



;a*a 



Ho = UHoU+ =Ho + vpS'^iO) (21) 

H,„t ^ UH,U+ = -§5-^(0) + -|-5-cos$,/(0) 

Hence, for the special value of the couphng, J^ — 
2'kvf , ^s no longer enters the hamiltonian, which 
takes the very simple form: 



Ho + —S'=cos^sf{0) 

TTflo 



(22) 



This describes a resonant level model for the 
Majorana fermion x^ = cos^sf, with scattering 
at the Fermi level. So x^ (which was left un- 
changed by the rotation U) suffers a tt/2 phase 
shift or, more appropriately for a real object, suf- 
fers a change of boundary condition from periodic 
{X''{0+) = X^(0")) to antiperiodic (x''(0+) = 
~X^(0 ))• The 'twist' operator cr^ which achieves 
this change of boundary conditions when acting 
on x^ ha-s dimension A{az) = 1/16, correspond- 
ing to the shift in the ground-state energy un- 
der a change of boundary conditions of one free 
Majorana fermion: LAE/irvp = 1/16. Hence 
A(a,) = 1/16. 

In contrast, the transformed fields x^ a-^d x^ 
drop out from the hamiltonian at the solvable 



point, and hence remain free and unaffected. The 
relation of these transformed fields to the original 
ones is easily worked out by considering the Dirac 
fermion: 

^s ^ -T^e-^"' (23) 

and observing that: 



-iTTS'sign{x)jj r^\ 



(24) 



Since S^ has eigenvalues ±1/2, one sees from this 
equation that 'ips{0~^) = e^*'^Vs(0~) and hence 
that the original field ips also suffers a phase shift 
S = tt/2 at the solvable point. In this case, the 
operator U which achieves this has the simple bo- 
son representation (g^ because it changes both 
components of a Dirac fermion. Its dimension is 
thus 1/8 = (5^/27r^, twice that of a twist operator 
for a single Majorana fermion. 

Hence, the crucial feature of the Emery- 
Kivelson solvable point is that the three (original) 
Majorana fermions ;^^'i''^ introduced in Eq.([15|) in 
order to represent the total spin current suffer a 
change of boundary condition from periodic to an- 
tiperiodic. (Actually, as discussed below, we need 
to consider both the antipcriodic and periodic 
Majorana sectors to construct the free fermion 
spectrum of the four original Dirac fermions with 
e.g antipcriodic boundary conditions. At the 
Kondo fixed point, the Majoranas which were 
originally periodic become antipcriodic, and the 
Majoranas which were originally antipcriodic for- 
mally become periodic, but in fact do not con- 
tribute to the spectrum at the Kondo fixed point.) 
This interpretation of the boundary condition as- 
sociated with the non-trivial fixed point has also 
been discussed independently by Maldacena and 
Ludwig Ig] and more recently by Ye g. In 
the bosonic language, the new boundary condi- 



tion reads: ^^(O^ 



$,(0-) , $,/(0^ 



IT — $5/(0"). The operator which achieves this 
change of boundary conditions has been con- 
structed in the present formulation (which is not 
explicitly SU{2) symmetric) as the product az -U 
of the transformation hi and of the twist oper- 
ator for x^- The total dimension of this opera- 
tor is thus 1/16 + 1/8 = 3/16. Within the CFT 



approach (which is manifestly SU{2) invariant) 
this operator is identified as the primary opera- 
tor corresponding to the fundamental representa- 
tion j = 1/2 of the SU2{2)^ Kac-Moody algebra, 
i.e to the WZW field ga- Note that its dimen- 
sion is indeed j{j + l)/4 = 3/16. Thus, the ex- 
plicit solution using abelian bosonisation at the 
Emery-Kivelson solvable point provides us with 
an explicit derivation of the 'fusion principle' of 
the CFT solution. 

1.4. Finite-size spectrum 

The whole finite-size spectrum can be simply 
recovered from the Majorana fermion approach 
(see also [||). As in the CFT approach, the first 
task is to construct the spectrum at the free- 
fcrmion fixed point. Let us choose antipcriodic 
boundary conditions for the four original Dirac 
fermions ipia- The corresponding free-fermion 
partition function reads: 



jpF 



JJ(1 + 5^+1/2)8 



(25) 



where q = ^-^^^/l^'"!' _ -yyg nmst find how to con- 
struct this spectrum in terms of 8 free Majorana 
fermions, 3 of them (%) being associated with the 
5^/2(2) spin degrees of freedom, and 5 being asso- 
ciated with the 5*01(5) charge-flavour degrees of 
freedom. It turns out that both the antipcriodic 
and the periodic sector must be considered, to- 
gether with a projection onto states with an even 
total fermion number F = Fspin ~^Fc-f- The cor- 
responding partition function in then obtained in 
the form (using formulas of Ref . (t) ) : 



Z^^=tr^i±(_i)I,L, 






m\% 



q 



1/6 



I-2U/ n / '^ 2\\ n ) ^ 



)8} (26) 



in which the Jacobi functions 0's and ry have argu- 
ment q. This can be checked to coincide precisely 
with (H). 

In order to obtain the spectrum at the Kondo 
fixed point from this expression, one must per- 



form a modular transform q 



= p-T^Pvp/L 



and twist the boundary conditions of the three 
Majorana fermions associated with the spin sec- 
tor, as explained above. This amounts to act- 
ing on the states with the operator (— l)-^"?'". It 
turns out that only the sector which was origi- 
nally antipcriodic at the free-ferniion fixed-point 
(and which has been twisted to periodic) con- 
tributes to the spectrum at the Kondo fixed point. 
The resulting partition function (with ground- 
state energy substracted) reads: 



2 + 4u;i/^ + 10^1/2 



(v'tr(Vt)n 



"~ "" (28) 

12u;5/8 + 26w + - •(29) 



in which the 0's and 77 have argument w. 
From this, the lowest energy levels (and de- 
generacies) are read off as: LAE/irvp — 
0(2) , 1/8(4) , 1/2(10) , 5/8(12) , 1(26), • • •. This 
can be checked to coincide with the result ob- 
tained from the fusion principle (0) in the CFT 
approach. 

1.5. Physical properties at and near the 
fixed point. 

Physical properties associated with the stable 
fixed point can also be analyzed in the bosonised 
language, and connections made with the CFT 
approach. It is convenient 0| to fully refermionize 
the hamiltonian in terms of the Dirac fermion ips 
and its spin- flavour counterpart: 



i'sf^ 



1 



-i<E>s f 



\/2TTao 



(30) 



and to represent the impurity spin components in 
terms of three local Majorana fermions: 



g^ 



V2' 



Sv 



V2' 



S"" 



-iab 



(31) 



The hamiltonian can be written as a sum of four 
terms H = He + Hf + Hgf + Hs , where He and 
Hf are the free charge and flavour parts, and the 
spin-flavour and spin parts read, in terms of these 
fermionic variables: 



H. 






(32) 



Hs = ivF f^^ dx'ipl{x)^'ilJs{x) 



(33) 

where the coupling constant X = Jz — 2'KVp mea- 
sures the deviation from the solvable point. The 
interpretation of Hgf as a resonant level model is 
very clear in this form. One also sees that only the 
b component of the impurity spin degrees of free- 
dom hybridizes at the solvable point, in the same 
way that only the i/jL- -|- ipsf oc Xz component of 
the conduction electron enters Hsf. a remains a 
local fermion, with non-decaying correlations: 



<Ta(0)a(T)>=-isign(T) 



(34) 



This gives rise to a residual entropy at the fixed 



point Si, 



ln2/2 



In a sense, only half of 



the impurity and conduction electron spin-flavour 
degrees of freedom hybridise. In contrast to a, b 
hybridises with x^ ^-nd thus acquires dimension 
1/2, so that impurity spin correlations at the fixed 
point decay as: 

< S{0)S{t) >cx< 6(0)6(r) >cx l/r (35) 

Exactly at the solvable point (where one has a 
quadratic fermion hamiltonian) , the impurity sus- 
ceptibility Ximp vanishes and the specific heat co- 
efficient C/T is non-singular H. The non- Fermi 
liquid singularities in these quantities are con- 
trolled by the first non-trivial order of perturba- 
tion theory in the coupling A away from the solv- 
able point. Indeed, the corresponding operator 
a6^+(0)?/'s(0), of dimension 3/2, precisely coin- 
cides with the CFT identification of the leading 
irrelevant operator at the non-trivial fixed point. 
To see this, we first note that b can be traded for 
X^ in the expression of this operator since the two 
fields hybridize at the solvable point. Then, we 
note that 4't4's oc X^X^ 1 so that the operator can 
be written as: 

i<^abcX''x''x'' oc (x X x).x (36) 

The triplet of Majorana fermions Xj of dimension 
1/2, is the (j>j=i primary operator of the Kac- 
Moody algebra S'[/2(2), so that (p9) appears as 
the descendant operator J7_i ■ 4>j=ij product of 
this operator and the spin current, which is pre- 
cisely the CFT identification. Perturbation the- 
ory in this operator |§,@i|| leads to the non-Fermi 



liquid behaviour C /T ex Ximp oc InT with a uni- 
versal Wilson ration Rw = 8/3. 

It is important to note that the original impu- 
rity spin variables (before the transformation lA 
is performed) also have a simple representation 
in terms of the Majorana ferniions x ^-nd a, h. In- 
deed: 

5+ EE 5^ + iSy = U+S+U = e''^''*'-(") [b + i a)e~''^''*' 
= (6-fia)V^,(0)/V2 

and 5^ = S*^. Hence, we obtain: 

S-^ = ax^'lO) - &X^(0) , Sy = ax^(O) + fox'^W 

S'^=ia6 (37) 

The most relevant part of the impurity spin op- 
erator, of dimension 1/2 can thus be identified 
with a x(0) (where the fermion b has been traded 
for x^lO)? with which it hybridizes). This identi- 
fication plays a crucial role in the solution of the 
two-impurity problem presented below. 

2. The two-impurity, tv^ro-channel Kondo 
model 

We now turn to the 'Kondo quartet' model, of 
two spin-1/2 impurities coupled to two channels 
of conduction electrons. This model was first in- 
vestigated using numerical rcnormalization group 
methods |10|-p^. Recently, an analytical solu- 
tion of the low-energy universal properties of the 
model was found (@ and ||l|-|lll). 

In its original formulation, the problem involves 
two spin-1/2 impurities, Si, S2 located at posi- 
tions ±i?/2 and interacting with two-channel of 
free conduction electrons. The electron gas is de- 
scribed by operators Cj:^^ , where k is the momen- 
tum, i = 1, 2 is the channel index and a =1 , 1 the 
spin index. The hamiltonian reads: 



H = Efc„e(fc)4^c,~„ + MSi-H+m) + 



S2-si-R/2)]-ISi-S2 



(38) 



In these expressions, s{r) is the total conduc- 
tion electron spin density at position r, Jk is the 
Kondo coupling (which we always take to be an- 
tiferromagnetic Jk > 0) and the inter-impurity 
coupling / has been taken as an independent pa- 
rameter. Of course, such a coupling is generated 



anyhow from the second-order process involving 
the conduction electron bath {i.e the RKKY in- 
teraction, of order J]^ for small Jk)- Standard 
renormalization-group arguments allow us to in- 
clude this coupling by hand (see also |1^ ). 

The problem ( p8|) can be considerably simpli- 
fied by reducing it exactly (at low energy) to 
r^ one-dimensional model involving 8 left-moving 
(chiral) fermionic fields ipiia{x), where the extra- 
index I = 1,2 originates from the presence of the 
two impurity sites. The procedure has been de- 
scribed in detail by other authors (see 6.17 [|9|), 
and we shall only quote the final form of the 
hamiltonian here: 

H = ivF J2ua J-^r dxil^ljx)£'iPua{x) 

+j+{Si + s2).iJiio)+j2m 

-HJ™(^l-^2).(ji(0)-^2(0)) 
+ J_(5i + 52).E,,a/3(^L(0)¥V'2.M0) 

+V'L(o)¥v'H/3(o)) 

-IS1.S2 (39) 

In this formula Ji{x) = Y.i,a(i^\ia(^)^'^iip{^) 
denotes the spin-current at position x of conduc- 
tion electrons corresponding to I = 1,2. The 
reader is directed to Ref.fl^ for the expression 
of the (bare values of) the coupling constants 
J+, Jm, J- in terms of Jk and R. Alternatively, 
one could work in the even-odd basis: 

Ipeia = (V'l,ia+V'2,m)/V2 , Tpoia = ('/'l,ia-'02,ia)/ ^2(40) 

in terms of which the Kondo couplings in Eq. (|39| ) 
read: 

{Si + S2)-{rJM + ^oJom 

+r™ (5i - 5^2) . E.,a/3(^L(0)¥V'o.M0) + h.c) 

Where the couplings r^o.m are the ones used in 
Ref.|0|: 

Te = (J++J-)/2 , r„ = (J+-J-)/2 , r,„ = J™(4i) 

The hamiltonian (p9|) is invariant under a par- 
ity transformation, which exchanges the indices 
1 — 1,2 for both impurity spins and conduction 
electrons: 



Parity : ipua ^ ip2ia , 5*1 <-» 5*2 



(42) 



Exchanging the even and odd combinations: 



■0e 



^o,ja ; ^1 



V'l.ja , 4'2,ia ^ -■02,ia(43) 



is a discrete symmetry of all terms in (p^ , except 
the J_ oc Fe — To coupling. 

When written in the simplified form (p9|), the 
hamiltonian is automatically invariant under a 
particle-hole transformation. Breaking particle- 
hole symmetry in the original hamiltonian of the 
model will induce potential scattering terms at 
low-energy, of the form: 



Vei'ti'e + Vo^ + i>o 



(44) 



2.1. Overview of the results 

We give here a brief overview of the results 
that will be derived in the following. Let us 
start by discussing the possible fixed points of the 
model in the space of the various coupling con- 
stants. There are several couplings {I,J±,Jm) 
in the hamiltonian (g^) , but we shall show below 
that the whole phase diagram can actually be dis- 
cussed in terms of only two parameters. These 
parameters (denoted x and y in the following) 
are phase-shifts which depend in a non-universal 
way on the original coupling constants, and can 
only be calculated analytically in some limiting 
regime. These two phase shifts are the natural 
variables in terms of which the problem is best 
described, and in terms of which numerical re- 
sults should be interpreted. The precise phase 
diagram in terms of these two parameters will be 
found in analytical form in Sec. 2.3. 



Here we follow a route which is physically more 
transparent, and display a schematic version of 
the phase diagram in the space of two variables, 
in terms of which the phase diagram has been 
investigated numerically by Ingersent and Jones 
11,1^]. These two parameters are the ratio I/Tk 
of the RKKY interaction to the Kondo temper- 
ature of the single impurity model and the ratio 
(Fe — Fo)/((Fe + Fo) ex J_ measuring the rela- 
tive asymmetry between the odd and even Kondo 
couplings. In order to reduce oneself to this pa- 
rameter space, one possibility is to assume that 
the coupling 5J ^ J+ — Jm is taken to be zero 
(we shall see that 6 J actually corresponds to an 
irrelevant operator around the decoupled impu- 
rity fixed point). 



In this restricted two-parameter space, the 
schematic phase diagram of the model is depicted 
in Fig.|l|. The figure is not meant to be quanti- 
tatively precise (a numerical determination using 
the Wilson RG has been given by Ingersent and 
Jones in [|ri|,Q), but merely indicates the var- 
ious possible fixed points and directions of the 
RG flow. Some of the fixed points in Fig.^ cor- 






X 



/ J . 



y^^lT 









\ 



^ 



X 






\ 



\ 



I/T, 



Figure 1. Schematic phase diagram in the 
[I/TK.J-/J+ = (Fe-Fo)/(Fe + Fo)) plauc. The 
marginal domain (made up of a surface and half- 
line) is indicated in dark. Outside this domain, 
the fiow is towards Fermi-liquid fixed points cor- 
responding either to even (or odd) Kondo effect 
or to the free-fermion (=RKKY) fixed point. The 
phase diagram is symmetric under J_ -^ — J_. 



respond to simple limits which can be easily dis- 
cussed qualitatively ||lO|] : 

- The fixed point at / = 5J = 0, J_ = 
(Fe = Fo) corresponds to two decoupled impu- 
rities interacting with two channels of conduction 
electrons. The Nozieres-Blandin single impurity 
physics is found at this fixed point. 

- The fixed point at / = -oo, J_ = (Fg == F^) 
obviously corresponds to free electrons. Indeed, 
for very strong antiferromagnetic RKKY cou- 
pling, the impurity spins bind into a singlet state, 
leaving a local Fermi liquid with no residual scat- 



tering for both the even and the odd combinations 

(<5e = <5o = 0). 

- The fixed point at / = +00, ,/_ = (Fe = To) 
corresponds to binding the impurity spins into 
an S" = 1 triplet state, while maintaining per- 
fect symmetry between the odd and even parity 
combinations. As conjectured in [[lO[ , it is natu- 
ral to expect that the physics at this fixed point 
is identical to that of a single S — 1 impurity 
with twice two channels of conduction electrons, 
i.e four channels. Since 4 > 2x1, this is an 
overscreened non- Fermi liquid fixed point, with 
II : X^mp '-' C/T -- T-i/3. We shaU prove this 
conjecture in the following. 

- Finally, the fixed point at / = 0, J_ — > +00 
(Fe >> Fo) has also a simple physical interpre- 
tation. At this fixed point, the odd electrons de- 
couple and the Kondo effect takes place for each 
impurity with the even combination only. Hence, 
the scattering is characterized by Se = 7r/2, So = 
0). Again, one has a local Fermi- liquid. When 
the RKKY coupling is ferromagnetic (/ > 0), we 
find that an arbitrary small odd-even asymme- 
try drives the system towards this fixed point. 
In view of the above interpretation of the strong 
ferromagnetic fixed point, we see that this even 
Kondo fixed point must have identical physics to 
a single-impurity, spin-1, 2-channel fixed point, 
which is indeed an exactly screened situation and 
hence a Fermi-liquid. On the antiferromagnetic 
side, this fixed-point controls only part of the pa- 
rameter space, namely large enough J_ and small 
enough / (Fig.|l|). 

Apart from these fixed points that can be found 
from simple arguments, the model also displays 
a continuous family of non-Fermi liquid fixed 
points (in dark on FigJ^). As indicated on Fig.n^, 
this continuous family is made up of a marginal 
line restricted to the J_ = (Fe = Fo) axis 
when / > is antiferromagnetic, and of a two- 
dimensional domain when / < is antiferro- 
magnetic. This results from the finding that the 
RKKY coupling is an exactly marginal perturba- 
tion around the decoupled impurity fixed point 
/ = J_ = 0, while the even-odd asymmetry is 
relevant for / > 0, and again exactly marginal for 
/ < 0. Within this marginal domain, universal 
low-energy properties {i.e the low-temperature 



critical behaviour of the various physical quan- 
tities and the form of the finite-size spectrum) 
depend continuously on two parameters, accord- 
ing to formulas derived below. As an example, 
for ferromagnetic RKKY and J_ = 0, a power- 
law behaviour Ximp ^ C/T ^ \/T^ is found, with 
an exponent varying continuously with I /Tk be- 
tween the limits 9{I = 0) = (corresponding to 
the Nozieres-Blandin logarithmic behaviour) and 
9 = 1/3 (corresponding to the S = 1, 4-channel 
model) . A separatrix exists in the phase-diagram, 
which separates the attraction basin of the free- 
electron (RKKY) fixed point and that of the even- 
parity Kondo fixed point. This separatrix hits the 
boundaries of the marginal domain at a multicrit- 
ical point (Figj^). 

Hence, the main physical message is that inter- 
impurity effects do not destroy non-Fermi liquid 
behaviour in this model (in the particle-hole sym- 
metric case), but do modify it as compared to 
the single-impurity case and generate a rich vari- 
ety of new behaviour. The competition between 
the Kondo screening and the RKKY ordering 
is highly involved in this model. This is to be 
contrasted with the two-impurity single-channel 
model (again at particle-hole symmetry). In this 
case, a RKKY (free electrons) fixed point con- 
trols the physics for sufficiently antiferromagnetic 
I /Tk < (I/Tk)c, while a phase with separate 
Kondo screening of the impurities is found for 
I/Tk > {I/Tk)c- Both of these fixed point are 
Fermi liquids and only the (unstable) multicrit- 
ical point at {I /Tk)c displays non-Fcrmi liqiud 
properties. 

2.2. The odd-even symmetric case 
2.2.1. Symmetries 

We concentrate in this section on the case 
J_ = (Fe = Fo) corresponding to a hamilto- 
nian invariant under odd-even exchange. When 
the coupling between impurities is turned off (/ = 
0, J-i- = Jm), the model has a global invariance 
{SU{2)spin ® SO{'b)fcY. At the decoupled im- 
purity fixed point, the symmetry algebra consists 
in two copies of a product of Kac-Moody alge- 
bra for spin, channel and charge: (S'f/2(2)s (g) 
SOi[^)fcY When coupling the two impurities 
(/ 7^ 0, J+ 7^ Jm) while keeping J_ = 0, the 



independent charge and flavour symmetries are 
left unaffected (this is no longer true for J_ 7^ 0), 
but the spin symmetry is reduced to the 'diago- 
nal' SU{2) for which ^1 and '4'2 are transformed 
identically. This reduces the conformal symme- 
try of the spin sector at a fixed point to a SU4{2) 
algebra. The generators of this algebra are the 
sum of the generators of the two SU2{2)s for 
each impurity, that is the sum of the spin currents 
Jiix) -\- J^ix). Hence we must understand how 
the product SUi(2)s ® 5'C/2(2)s can be decom- 
posed into 5[/4(2)s plus some residual degrees of 
freedom. This procedure is known in CFT as a 
'coset construction' |20|. The residual degrees of 
freedom define an algebra A(2, 2) such that: 

5C/2(2), ® 5772(2), =51/4(2), (»A(2, 2) (45) 

The conservation of the total number of degrees of 
freedom between the two sides of (Esl) is expressed 
by the conservation of the total central charge. 
Since the central charge of 5[/2(2) is c = 3/2 
and that of 5^/4(2) is c = 2, we deduce that the 
central charge of the A(2, 2) algebra is c = 3/2 -|- 
3/2 — 2=1. This algebra is actually known from 
CFT [M,pl| to be a A^ = 1 superconformal unitary 
model corresponding to the ?7i = 4 member of 
the discrete series with central charges c = |(1 — 



m(Tn-f 2) 



2.2.2. Bosonisation 

We shall give a schematic derivation of the 
operator content of the algebra A(2, 2) using 
the Majorana fermion representation and abelian 
bosonisation. We introduce two triplets of Majo- 
rana fermions Xi'2'^ such that the 5C/2(2) spin 
currents Jiix) {I = l,2;a = x,y,z) read: 
Jf'ix) = ieabcx\xi- These six Majorana fermions 
can be combined into three Dirac fermions, and 
bosonised as: 



Xl{x)+ix2{x) 



1 



^/TkOq 



'*"(^) ; a = x,y,z{AQ) 



In terms of these fields, the total spin current cor- 
responding to the diagonal SUi{2) algebra reads: 

J-~Jf + J^ = cos{^y - $,) , jy = cos{<^^ 



It is convenient to introduce three linear combi- 
nations of boson fields as follows: 

^ =-!=($, + $^-2$,) (47) 

Let us also note for further use the inverse rela- 
tions: 

In term of these combinations, the components of 
the total spin current read: 

J- = cos(ii^) , jy = cos{ii+^) 

J"' = cos(^/2^) (49) 

Note that $ does not enter these expressions. 
The two bosons /.t, v are associated with the 
51/4(2) algebra (which has c = 2), while $ cor- 
responds to the residual A(2, 2) degree of free- 
dom (c = 1). This boson is compact, with a ra- 
dius R — \/3, which means that $ and <& -I- 27ri? 
are identified. Hence, A(2, 2) contains all opera- 
tors of the form e±^("v^+™/2N/3)* ^nd 9"$ (with 
n, 771 integers). This yields primary operators 
of dimensions: (0) , (1/24) , (1/6) , (3/8) , (1/6 + 
1/2) , (1), as summarized in Table ^. In addition, 
it can be shown that $ and — $ must be identi- 
fied, so that ^(2, 2) is in fact an orbifold theory. 
This implies that the algebra contains, in addi- 
tion to the operators above, two twist operators 
of dimension 1/16 and two twist operators of di- 
mension 9/16. The full set of primary operators 
of the A{2, 2) algebra and their boson represen- 
tation (when it exists) is given in Table ||. 

At this stage, we can compare the coset con- 
struction made above with that encountered in 
the 2-impurity, single channel Kondo model [|l8| . 
There, the symmetry algebra at the decoupled 
impurity fixed point is SUi{2) (g) SUi{2). It is 
broken down to SU2(2) when coupling the two 
impurities. The coset algebra for the residual de- 
grees of freedom thus have central charge c = 
1 + 1 — 3/2 = 1/2 and identifies with the Ising 
<j) ynodel: 

SUi{2) ® SUi{2) = SIM2) ® (Ising) (50) 



W 



(48) 
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Table 1 

Primary operators of the A{2, 2) algebra. (A) labels an operator of dimension A. The second line displays 
the bosonic representation of the operator (when it exists). NS and R stand for the Ncveu-Schwarz and 
Ramond sectors of the algebra. 



(0) 


(^) 


ij-e) 


i^e) 


a) 


(1) 


i^e + k) 


m 


a + k) 


(1) 


1 


e^'^ 






e±^* 








e±^« 


a$ 


NS 


R 


NS 


R 


NS 


R 


NS 


R 


NS 


NS 



The Ising algebra contains three primary opera- 
tors, of dimension (identity), 1/16 (order pa- 
rameter), and 1/2 (energy density). Affleck and 
Ludwig lia] have shown that the properties of 
the non-Fermi liquid unstable critical point pTJ I 
encountered for this model are obtained from the 
decoupled impurity fixed point, by a fusion in- 
volving the order-parameter operator of the Ising 
sector. In the present case, we shall see that no 
direct fusion principle with an operator of the 
A{2, 2) coset exists. We shall nevertheless de- 
rive a generalization of this principle involving 
boundary-condition changing operators which are 
not in general primary operators of the coset al- 
gebra, except at special points in the phase dia- 
gram. 

2.2.3. The RKKY interaction as an exactly 
marginal perturbation 

We consider the effect of turning on the cou- 
plings I and 5J = J+ — Jm away from the decou- 
pled impurity fixed point (but keeping J_ =0). 
At the decoupled point, the impurity spin has di- 
mension 1/2, and the conduction electron spin 
current has dimension 1. Thus, the RKKY cou- 
pling has dimension 1 and is marginal at lowest 
order, while 5J has dimension 3/2 and is irrel- 
evant. We shall see that / is actually exactly 
marginal to all orders, while the dimension of 5 J 
is actually continuously modified as one departs 
from the decoupled fixed point. In order to rep- 
resent the RKKY term, we make use of the two 
sets of Majorana fermions xi, X2 above and intro- 
duce, as in Sec|^, two sets of Majorana fermions 
O'Ubi {I = 1,2) in order to represent the impu- 
rity spins Si as in Eq.nq. As explained at the 



end of section |l|, the operator of lowest dimen- 
sion {— 1/2) contained in Si reads aixi in terms 
of these variables. Hence, the perturbing term of 
lower dimension (= 1) associated with the RKKY 
interaction reads: 



dtaia2Xi-X2 



(51) 



In the bosonic language above, this translates into 
an induced boundary term in the A(2, 2) sector 
of the hamiltonian: 

where we have set: 

d+ = (ai -f ia2)/V2 (53) 

with a normalisation a^ — a\ — 1/2, so that 
{d,d^} = 1. Here, / is some (non-universal) 
function of / and dJ=J+ — Jm, which has a per- 
turbative expansion I — I + 0{I^,5,P) around 
the decoupled point. Note that the perturbing 
term in Eq.(p3) is perfectly compatible with the 
orbifold nature of the A{2, 2) algebra, provided 
the transformation $ — > — <I> is always made si- 
multaneously with d — > d"*" on the local degree 
of freedom. Since d remains a local fermion 
(with non-decaying correlations), it is clear from 
Eq.(|5^) that the RKKY coupling is associated 
with a dimension 1 operator and is thus an ex- 
actly marginal perturbation. This results in the 
line of fixed points on the Fg — To axis of Fig.|l|. 
The hamiltonian (|5^) is very similar to the 
X-ray edge problem in the bosonised form |2^. 
The interaction term on the boundary can be ab- 
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sorbed into a redefinition of the field: 

1, 



^{x) -^ $(a;) 



I—id+d 

Vf 



sign{x) 



(54) 



This changes the boundary condition on the field 
$ by a phase-shift 6, such that <1>(0+) - $(0~) = 
±26 when d^d — 1/2 = ±1, respectively. The 
weak-coupling expansion of this phase shift is: 



S 



I 



I 



2vp 2vp 



+ OiI^,6J^) 



(55) 



As observed by Schotte and Schotte p2| and 
also recently discussed in [p3|, this means that 
the interacting hamiltonian is related to the non- 
interacting one by a canonical transformation: 



H 



A{2,2) 



UtHj^^Us , Us 



„»f (d+d-i)<I.(0) 



(56) 



This transformation changes the dimension of the 
operators e*^* from ^ to ^(g ± |)2. Note that 
these dimensions do not correspond in general to 
dimensions of existing operators in the A{2,2) 
algebra. This reflects the fact that Us is itself 
not a primary operator of this algebra in general. 
Since the available dimensions in this algebra are 
given hy g = nv3 -I- m/2v3, the maximal possi- 
ble scattering (corresponding to a shift of unity 
of the integer m) is reached for ""^ = -^. 
Each fixed-point along the marginal line in fig.l3 
corresponds to a specific value of the phase-shift 
—6 max < (5 < 6 max-, which spccifics Completely all 
universal properties at the fixed point. It is con- 
venient to make use of the normalised parameter: 



2^/?>8 



(57) 



Note that x > (resp. a; < 0) corresponds to 
a ferromagnetic (resp. antiferromagnetic) RKKY 
coupling. As will be shown below, the marginal 
lines terminates at infinitely strong RKKY cou- 
pling / — +00 on the ferromagnetic side, corre- 
sponding to maximal scattering a; = -1-1, while 
the marginal behaviour only persists for x > x = 
Xmin = ~(3 — \/6) — —0.55 on the antiferromag- 
netic side. 

2.2.4. Physical properties 

We can investigate the operator content and 
low-temperature critical behaviour of physical 



quantities for non-zero /, bJ (corresponding to 
a specific value of a;), using the fact that a given 
operator O is changed to UsOUg' for a non-zero 
X. We shall deal first with impurity spin corre- 
lation functions. As already noted, the most rel- 
evant operator (of dimension 1/2) corresponding 
to each impurity spin at the decoupled impurity 
fixed point reads Si is aixi, and we thus have the 
identification: 



Sf-S^(xde-"^-+d-e 



h d+e-'^" 

**» , a^x,y,z (58) 



We act on these operators with the transforma- 
tion Us, using UsdUg' = exp (— ia;$(0)/-\/3) d. 
Hence the dimensions A^'" of Si ± 52 read, for 
arbitrary x: 



A Spt7l 



r+-. 



^spi'. 



1 , 1 + a: 



)'+3(59) 



From this, we find the low-temperature behaviour 
of the impurity spin susceptibility: 



\imp ^ j 



1/T 



dt 



TOM 



0{x) 



1-2A 



spin x{2—x) 



{x>{)) 



(60) 



while Ximp remains finite on the antiferromag- 
netic side X < 0. Note that the critical behaviour 
( |6C| ) interpolates continuously between the known 
limits ^: Ximp ^ ln{T) at the decoupled impu- 
rity fixed point, and Ximp ^ 1/T^'^ of the 5 = 1, 
four channel model at the strong coupling fer- 
romagnetic fixed point. Similarly, we find that 
the staggered susceptibility Xst (defined as the 
response to a field coupled to 5i — 52) is finite on 
the ferromagnetic side a; > 0, and diverges on the 
antiferromagnetic side as: 



Xst 



T0(- 



e{-x) = 1-2A 



Spt7l 



\x\i2 + x) 



(61) 



This continous dependance of the critical expo- 
nents on X establish the existence of a line of fixed 
points extending on both sides of the decoupled 
impurity fixed point. In order to find the precise 
extension of this line, and the low-temperature 
behaviour of the specific heat along it, we look 
for the leading irrelevant perturbation(s) compat- 
ible with parity and the even/odd symmetry. At 
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the decoupled impurity fixed point, there are two 
such operators, of dimension 3/2, corresponding 
to the couphngs J+ and Jm'- 



o+ = {Si + s2).{jiio)+j2m 

Ora={Sl-S2).{Jl{0)~j2m 



(62) 



We can find the bosonised form of 0+, using the 
representation of 5*1 + 5*2 given above , and the 
expression (B9) of the total spin current: 

Acting with Us on this expression, we find that 
both terms give rise for arbitrary x to an operator 
of dimension: 



^-'h^Z-w''''^"''' 



(64) 



Similarly, for Cm, we used the bosonised expres- 
sions: J7f — J2 ^ cos{^y + $^) (and permuta- 
tions) to obtain: 



Om oc d e 



-.^* + de^^Tf*--^^-;*'') + h.c (65) 



Upon acting with Us, the first operator becomes 
of dimension (\/3 + a;/V3)^/2, while the second 
has the same dimension than 0+. Hence the op- 
erator of lowest dimension which is generated has: 



A„, =mm{i(V3 + ^)2, A+} 



(66) 



The low-temperature behaviour of the specific 
heat is related to the dimension A of the irrel- 
evant operator of lowest dimension by: 



Hence: 

C/T ^T-'^ ,a = 3-2A 



(67) 



(68) 



For X > (ferromagnetic side), 0+ has the low- 
est dimension, and we find a = 3 — 2A-|_ = 
l„2A!_P™ = 0(x),sothat: 



X > : C/T ^ x»mp ^ r-^(") 



(69) 



and we expect a universal (x-dependent) Wilson 
ratio, which can be calculated using perturbation 



theory along the lines of Ref . [^ . For a; < how- 
ever (antiferromagnetic side), Om has the lowest 
dimension, and we find: 



<0: C/T - T-"(^) , a(a:) = 3-2A„ = 



,& + X 



(70) 



Note that Xst behaves with a different exponent. 
The limiting value of a; = Xmin, associated with 
the termination of the marginal line for antiferro- 
magnetic coupling, is reached when Om becomes 
relevant, i.e A„j = 1. This yields: 



Xmin = -(3- \/6) 



-0.55 



(71) 



For X < Xmin, the system flows to the 'strong 
antiferromagnetic' (or 'RKKY') fixed point, and 
the two impurities bind into a singlet state. Obvi- 
ously, this fixed point has the properties of a local 
Fcrmi-liquid. On the ferromagnetic side however, 
the marginal line extends all the way up to the 
infinitely strong RKKY coupling / = -|-cx) (corre- 
sponding to maximal scattering x = -1-1), since 
O-^ is still irrelevant at this point. Note that 
the critical properties derived above coincide, for 
X = +1, with that of the spin-1, 4-channel Kondo 
problem (with spin dimension 1/3 and leading ir- 
relevant operator of dimension 4/3), in agreement 
with the conjecture made in Ref.||in] and with the 
physical picture that the two impurities bind into 
a S* = 1 triplet state with effectively 4 channel of 
conduction electrons (because of even-odd sym- 
metry) . 

2.2.5. Finite-size spectrum 

We now investigate the finite-spectrum of the 
model at a given fixed point along the marginal 
line, as a function of x. In order to derive the 
spectrum, one first classifies the states of the 
decoupled impurity fixed point according to the 
SiM2)^ (g) A(2, 2) (g) (50^(5)^^2 decomposition, 
and then act on each state with the transforma- 
tion Us- This modifies the contribution of the 
A{2, 2) sector to the total energy of the state. 
The dimension 1/16 of the twist operators can 
be shown to be unchanged by the action of Ug. 
Under multiplication by Us, the dimension of 
an operator e'*^* is changed to i(fc ± ^)^, for 
d^d — 1/2 = ±1/2, respectively. Hence, we also 
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need to associate to each state an eigenvalue of 
d+d — 1/2 == ±1/2 to decide which of the two 
possible new dimensions is produced. This can 
be done, when constructing the spectrum at the 
decoupled point, by keeping track of the relative 
sign between the impurity spin and the total spin 
of the state. More precisely, the impurity spin {S) 
is proportional to the adjoint operator of SU2{2) 
(x) up to a sign which depends on the state. The 
product of these two signs for ^ = 1, 2 yields the 
eigenvalue of 2d'^d ~ 1. In particular, this 'selec- 
tion rule' is essential to insure that the spectrum 
of the S* = 1 four-channel model is indeed ob- 
tained at a; = +1. This observation is also a way 
to recover this selection rule in an empirical man- 
ner. 

The resulting spectrum is displayed in Table 0. 
The normalised excitation energy of a given state, 
^^^ (with L the radial length of the bulk system 
and vp the Fermi velocity), is obtained from the 
total dimension Atot given in the table, by the 
formula: 

1 _ (1 - .X)2 

3 " 



LAE 



= At 



nvp 
LAE 



24 



TTVp 



= A.„,-^(l + |f 



(for x>0) (72) 
{for X <0) (73) 



The ground-state is the triplet of lowest energy for 
X > (ferromagnetic coupling) , and the singlet of 
lowest energy for a; < (antiferromagnetic cou- 
pling) . Comparison with the numerical renormal- 
ization group results of K.Ingersent and B.Jones 
jlffl-E2| reveals that the numerical spectra are ex- 
cellently fitted by formula (uSh depending on the 
single parameter x (the agreement with the quan- 
tum numbers and degeneracies of each state has 
also been checked). Note that, despite the differ- 
ent degeneracies of the finite-size ground-state for 
X > and a; < 0, the residual entropy defined by 
taking the infinite-volume limit first is constant 
along the marginal line: 



Simv = lim lim ^S 



ln2 



(74) 



while Simp—0 at the strong antiferromagnetic 
fixed point. This is expected from the boundary 
version of the 'c-theorem'. 

The strong-coupling ferromagnetic fixed point 
found for a; = +1 deserves special comments. 



There, the impurity spins bind into an S* = 1 
triplet state. Thus, it should be possible to de- 
scribe this point, in the spirit of the CFT ap- 
proach |2], by a direct fusion in the spin SU4{2) 
sector with the triplet operator j — I (of di- 
mension j{j + l)/6 — 1/3) applied on the free 
fermion fixed point. One can check that per- 
forming this fusion gives precisely the spectrum 
of a single-impurity, four channel Kondo model. 
Furthermore, one can also check that exactly the 
same spectrum is obtained when acting on the 
spectrum of the decoupled impurity fixed point 
with the operator expz(d+(i — l/2)$(0)/'\/3 (cor- 
responding to Us for X = +1), taking into account 
the selection rule discussed above. Note that this 
operator is then a member of the A{2, 2) algebra, 
of dimension 1/24. Accordingly, the A{2,2) oper- 
ator labelling any eigenstate at x = is changed 
at a; = -1-1 into another operator of the A{2,2) 
algebra. This strong-coupling assignment is also 
displayed in Table g. It is compatible with the 
operator product expansion of the (^) operator 
with any other operator of the A{2, 2) algebra: 

(2l)x(2i)-(0) + (l) + (^) 

(s) X ije^NS -> {Te)R + (^)fl 

(54) X (ll)« ^ (ll)ws' + (11 + |)jvs 

(2l)x(l)-(i) + (2l) 

X (^)fl -^ (i^)a'S' + (]^ + h)NS 

' X (11 + ^)nS -> ( ]^)fl + (^)fl 

The selection rule discussed above dictates, for 
each state, which of the possible operators ap- 
pearing in the r.h.s of these rules are actually 
generated when going from the decoupled impu- 
rity spectrum to the strong ferromagnetic one. 
For example, the singlet state in Table and the 
multiplet of degeneracy 10 are both associated 
with the operator (|) for a: = 0, but are assigned 
different operators ((i -I- i) for the first one, (i) 
for the second) at the strong ferromagnetic point. 

2.3. The effect of odd-even asymmetry 

We now consider the effect of a non-zero value 
of the coupling J_, breaking the symmetry be- 
tween even and odd combinations (Fg ^ To). In 
the presence of this coupling, the independent 
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Tabic 2 

Finite-size spectrum of the first 76 low-lying states in the even-odd symmetric case J_ = (Fg = Fq). j is 
the total spin quantum number associated with S'C/2(4)^p™. The second column gives the charge-flavour 
representation in the 5*01(5)'^^ (E) SOi{5y^^ decomposition, a is the sign of d'^d — 1/2 = ±1/2. The 
third column displays the A{2, 2) operator associated with each eigenstate at the decoupled impurities 
fixed point (/ = 0, i.e a; = 0), whereas the fifth column displays the corresponding operator at the strong 
ferromagnetic fixed point (/ = -l-oo i.e x = +1). The degeneracy of each state is displayed in the last 
column, while Atot is the total conformal dimension at arbitrary x. The normalised excitation energies 
are given by LAE/ttvf = Atot - Af^t- 



j 


50(5)'=^® 50(5)'=" 


^(2,2) 
decoupled 


a 


Atot 


A(2,2) 
strong ferro. 


Deg. 





(1,1) 


(1) 


+ 


1(1 + f)^ 


(i + i) 


1 


1 


(1,1) 


(21) 


- 


1 , (i-^r 

3 ^ 24 
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flavour and charge global symmetry (50(5))" 
no longer holds: only identical transformations 
for both impurities are allowed. The symmetry 
group in the flavour and charge sector is thus bro- 
ken down to the 50(5) diagonal subgroup, the 
sum of the flavour and charge generators. Cor- 
respondingly, at a fixed point, the model has a 
5^/4(2)^ (g) 502(5),.^ Kac-Moody symmetry, and 
we must consider the following 'coset construc- 
tion' in the flavour-charge sector (in addition to 
the one above in the spin sector): 



50i(5)®50i(5) == 502(5) ®X 



(75) 



It is crucial in this case to take into account the 
full 50(5) symmetry of the charge-flavour sector, 
and not only its SU{2)f (g) C/(l)c subgroup. In- 
deed, this excludes some marginal perturbations 
that would be naively admissible in the indepen- 
dent flavour/charge language and explains how 
multiplets of states at J_ = are broken into 
multiplets of smaller degeneracies when J_ 7^ 0. 
The current algebra of 50i(5) is essentially a 
straightforward generalization of that of SU2{2) 



to the case of five Majorana fermions (it has cen- 
tral charge c = 5/2). Again, five Dirac fermions 
are built out of the two sets of Majorana's, and 
converted into five bosonic fields <I>i. The diagonal 
combination <& = ^^ <I>i/-\/5 is associated with 
the coset algebra K, which has central charge 
c — 1. The four other combinations correspond 
to 502(5) (with c = 4). The coset X is a subalge- 

bra of the product of coscts A{2, 2)/ (g) C^(l)ci_c2' 
which have both c = 1. Accordingly, the bo- 
son $i is a specific linear combinations of the 
two bosons associated with the flavour and charge 
cosets, namely: 



$ 



$/ 



2 <^l 



$g 



V2 



(76) 



We now ask whether turning on a small cou- 
pling J_ is a relevant perturbation at a given fixed 
point on the marginal line of the previous section. 
To answer this, we consider the corresponding op- 
erator: 

0_ ^ (5i + 52).E.,„0(V'L(O)¥^2.MO) 
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+V'L(o)¥^H/3(o)) 



(77) 



2\ 2 



.)2 , 5 



A naive counting of dimension would suggest that 
this operator has dimension 1/2 + 1 = 3/2 at the 
decoupled impurity fixed point. This is incorrect 
however, since ©_ is a product of an impurity 
and fermion parts which are not independent. In- 
stead, we shall see that O- is, to lowest order in 
J_, a marginal operator (of dimension 1) at the 
decoupled impurity fixed point. It is convenient 
to find the bosonised form of 0_ at this fixed 
point. S1 + S2 is represented as in eg. (|58|) , and we 
need to concentrate on ■(/'iia(0)^^i/'2i/3(0) + h.c. 
This operator is a spin triplet, and a flavour 
and charge singlet. Using the spin- flavour-charge 
decomposition 5X72 (2)s «) 5772(2)/ (g) C/(l)c for 
each impurity index I — 1,2, we have to build 
a triplet combination out of two j — 1/2 opera- 
tors of the SU2{'2)s algebra (each of dimension 
j{j + l)/4 = 3/16, so that the spin contribu- 
tion to the overall dimension is 3/16 -I- 3/16 = 
3/8). Similarly, we build a singlet operator out 
of two j — 1/2 operators of the S'C/2(2)/ alge- 
bra (also of dimension 3/8), and a singlet charge 
combination (of dimension 1/8 + 1/8 — 1/4). 
If we represent each fermion i/j^l^^ using non- 



abelian bosonisation (Eq.(|6|)j, as Qii^ii-i 



this corresponds to 

9ta^c,0920h+h2^e^^^^-^'^r/2 



forming the combination 



The flavour singlet 



combination /ij'j/i2i is built out of the three oper- 
ators (of dimension 1/8) e'*"/^,a — x,y,z (which 
give a 7r/2 phase shift to one of the three flavour 
Dirac fermions), as expz(<i>/ + <&/ + ${)/2 = 
expi\/3$^/2. Similarly, we built the spin 
triplet giaffai3g2i3 out of the e**"/^ operators as: 
expi(— <&!, + ^y + <&z)/2) (and permutations), i.e 
as: expi(\/3$/2 — $a) ,a — x,y,z. Multiplying 
by Si + S2 given by Eq. (Bq) , we find that the op- 
erator of lowest dimension contained in O- reads 
in bosonised form: 



dexp{i(^$(0) 



+ - 



2 
f?(o)- 



^2(0) 



)} +h.c 



Hence, this operator is of dimension 1 at the de- 
coupled impurity fixed point, as announced. Un- 
der the action of Us , the spin part is transformed 
into expi{^/3/2 — x/^/3)^{0) and hence the di- 



mension is changed to: A^{x) — 2^-2 T?; 

Since A_ (x) < 1 for a; > 0, we conclude that the 
coupling J_ is a relevant perturbation on the fer- 
romagnetic side, and that the system flows away 
from the marginal line. For e.g Fg > Fq, the 
flow must be towards a fixed point where the odd 
combination eventually decouples, and thus cor- 
responds |1^ to an even-parity Kondo effect with 
S = 1 and two channel of conduction electrons, 
which is a Fermi liquid fixed point (cf. fig.??). 
On the antiferromagnctic side of the line (a:: < 0) 
however, wc find that A_(x) > 1 and thus that a 
small J_ is an irrelevant perturbation. Hence, the 
marginal behaviour must he preserved for small 
enough J_ on this side. 

In order to find the new marginal operators 
generated by the coupling J_ for x < 0, we con- 
sider the perturbation expansion in the operator 
d?^). This operator is irrelevant, but its oper- 
ator product expansion with itself will generate 
the marginal operator (d+d— l/2)9a;<I'(0) (induc- 
ing simply a renormalisation of /, proportional to 
J^ for small J_), but also the marginal operator 
(d+d — 1/2)92:'I'(0), where $ is precisely the lin- 
ear combination of flavour and charge bosons ( |7q ) 
corresponding to the coset K. This operator was 
forbidden for J_ — 0. The fact that precisely this 
combination appears is due to the 5*0(5) sym- 
metry of the theory. Had we used the separate 
flavour-charge decompositions of, we would have 
concluded incorrectly that two new marginal per- 
turbations dx^^{0) and dx{^i — ^2) were a priori 
allowed, corresponding to the two c = 1 cosets 

A{2,2)f and C^(l)ci_c2- Hence, we find that, in 
addition to the induced boundary term (52) in 
the spin sector, we have a single other boundary 
term induced in the flavour-charge sector, which 
reads: 



1 (9$ 



(79) 



The problem is again solved by a canonical trans- 
formation which now involves a single additional 
phase shift S' associated with (pa), depending 
on the couplings in a non-universal manner with 
6' oc J^ for small J_. In the following, we make 
use of a parameter y such that S' /tt = \/5y/2. 
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The overall canonical transformation reads: 



Uss> = exp{i((i^(i- 






$(0)- 



\/5y 



$(0))}(80) 



For antiferromagnetic RKKY coupling (a; < 0), 
the marginal line extends into a two-dimensional 
marginal domain parametrized by the two param- 
eters X and y. In order to find its precise exten- 
sion in the (x, y) plane, we discuss the operator 
content at a given fixed point of the marginal do- 
main, concentrating on the irrelevant operator(s) 
of lowest dimension. This analysis will also yield 
the low-temperature behaviour of physical quan- 
tities. The boundaries of the domain will be char- 
acterized by one of these operators becoming rel- 
evant. It turns out that we have to discuss only 
Om and 0-, all other operators (including 0+) 
being of higher dimension in this region of the 
phase diagram. We use the bosonised form ( JTq ) 
of O- to find the dimension of the transformed 
operator Uss'O-Ut^,, which reads: 



A_{x,y) 



6^2 



2^2 



y? 



Similarly, using eq.(p3) for Or, 



A™(x,y) = ^(l + |)2 



2^ 



(81) 



(82) 



Note that these expressions coincide with those 
established above in the limit y = 0. Thus, we 
find that these operators are irrelevant (A > 1) 
outside the ellipses of equations: A_(a;,y) = 1 
(for 0-) and A™(x,y) = 1 (for Om)- The lo- 
cation of these ellipses in the {x,y) plane is de- 
picted on fig.0. Remarkably, they are tangent 
at a single point {xc,yc) — (—3/4,1/4). The 
shaded area in fig.||, delimited by these two el- 
lipses and the two coordinate axis, together with 
the whole half-axis y = 0, a; > corresponds to 
the marginal domain. When crossing the bound- 
ary to the right of the point (xcj/c) ('ferromag- 
netic side') , O- becomes relevant and the system 
flows to the even-parity (or odd-parity) Fermi- 
liquid fixed point (except when y is strictly zero). 
When crossing the boundary to the left of {xc, yc) 
('antiferromagnetic side'), Om becomes relevant 
and the system flows to the strong antiferromag- 
netic Fermi-liquid fixed point / = — oo. Note that 




Figure 2. The marginal domain in the plane of 
the phase shift parameters (x, y) is made up of 
the shaded area and the half-axis a; > 0, y = 0. 
The (tangent) ellipses inside which the operators 
O- and Om become relevant are indicated: see 
text for a detailed explanation. 



y is even in J_ and thus that the phase diagram 
in the y > half-plane yields a symmetric phase 
diagram in the coupling J_, with odd and even 
coupling exchanged. Thus, we identify {xc,yc) 
as the multicritical point separating the attrac- 
tion basins of two different fixed points when the 
boundary of the marginal domain is crossed. It is 
worth noticing that at this point the dimension of 
the operator Uss' becomes x^/24-f5y^/8 = 1/16, 
hence suggesting that a direct interpretation of 
this point as a fusion from the decoupled fixed 
point with some twist operator is probably possi- 
ble. This is reminiscent of the multicritical point 
of the single-channel, two impurity model [|l8| . 

We also discuss the low-temperature behaviour 
of the impurity susceptibilities and specific heat. 
Acting on the bosonised form (58) of the opera- 
tors Si ± S2 with Uss' , we find their new dimen- 



A spin 



^spin _ I 
^- ^ 3 



6 
{1+xf 



5r 
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Hence, the impurity susceptibility is found to 
be finite within the whole antiferromagnetic part 
a; < of the marginal domain, while the staggered 
susceptibility (defined as above) diverges as: 

. |x|(2 + : 



Xst 



T- 



5y' 



(83) 



In order to find the behaviour of the specific heat 
inside the marginal domain, we have to find which 
of the two irrelevant operators Om and O^ has 
lowest dimension. Using Eqs.( |8l| , |8^ ), we find that 
A_ < Am to the right of the straight line 3a; + 
5y = 1 (also depicted as a dashed hue on fig.||), 
while Am < A_ to its left. Hence, the specific- 
heat behaves as: C/T ~ T"" with: 

a = 3-2A^{x,y) for3x + 5y<l,{yj^0) 
a = 3-2A™(x,j/) for 3x + 5y>l (84) 

The straight line 3x + by = 1 contains {xc,yc) 
(where A_ = A™ — 1) and crosses the y = 
axis at a; = —1/3 > Xmm- On this axis however, 
the operator O- is not allowed because of the 
odd-even symmetry, and the low-temperature be- 
haviour of C/T was found to be controlled by Om ■ 
Hence our results lead to the somewhat surprising 
conclusion that the exponent of the specific heat 
should change discontinuously for —1/3 < a: < 
when y is turned on. 

The finite-size spectrum of the model in the 
marginal surface can be obtained in a similar 
manner than for the even-odd symmetric case, by 
acting, for each state, on the operator in the coset 
A{2, 2)iS)K with the transformation Uss' ■ The re- 
sulting spectrum is also in good agreement with 
available numerical data, and will be reported in 
another article [mUi^. 

3. Conclusion 

We have shown that the competition between 
inter-impurity (RKKY) interactions and the (par- 
tial) screening associated with the Kondo effect, 
lead to a continuous two-parameter family of new 
non- Fermi liquid fixed points in the two-impurity, 
two-channel Kondo model at particle-hole sym- 
metry. An analytical solution of the universal 
low-energy properties of this model has been ob- 
tained, using a combination of bosonisation and 
conformal field theory techniques. 



Possible physical realizations of this model in- 
clude magnetic impurities in coupled Heisenberg 
spin chains |14| and coupled quantum dots de- 
vices |2jj2^. In those cases, one often faces an 
RKKY interaction which is not isotropic in spin 
space. This can be solved by bosonisation along 
the same lines as the isotropic case considered in 
this paper. The perturbation then acts not only 
in the coset algebra, but also in the SU 4(2) spin 
sector, and additional marginal operators (hence 
additional phase shifts) are induced. 

An outstanding question that has not been 
dealt with above is the effect of particle- hole sym- 
metry breaking on this family of fixed points. We 
have shown n& that it is a relevant perturbation 
near the strong ferromagnetic fixed-point (and 
hence, by continuity, along the whole marginal 
line for a; > 0). The situation on the antiferro- 
magnetic side is currently under investigation. 
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